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SHAPE DESIGN OF A FILLET 


Selection of the best shape of a fillet in a tension bar such that no 
yielding occurs has long attracted the attention of engineers. Dimensions and 
notations for the bar and fillet are shown in figure 1. Witlji symmetry, only the 
upper half of the bar is considered. The boundary segment T is to be varied, 
but with fixgd points at A and B. The segment T is the central line of the 
fillet and T and T are uniformly loaded edges. 

The optimal design problem is to find a boundary shape T ^ to minimize the 
total area of the fillet such that no yielding occurs. Constraints are placed 
on von Mises yield stress, averaged over small regions or finite elements ^ on 
which m^ is a characteristic function with value l/(area of and <f>(a(z)) is 
normalized von Mises yield stress. 

The classical boundary value problem is reduced to a variational or energy 
related problem which not only has excellent properties of existence and 
uniqueness but also provides the mathematical foundation for finite element 
analysis. The variational formulation may be viewed as the principle of virtual 
work and the finite element method as an application of the Galerkin method to 
the variational equation for approximate solution of the boundary value problem. 



Design Variable: Shape of T 1 

Cost: i|» 0 = J/ Q d£2 

Constraint: \|^ = //^ <()(a(z))m^ dft < 0 , k = 1,2,...,NE 

Virtual Work Equation: 

2 2 
a(z,z) = // l [o 1J (z)e J (z)]dfl = / 2 1 T z dT , 

i,j=l i=l 

for all kinematically admissible virtual displacements z 

Figure 1 
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MATERIAL DERIVATIVE AND ADJOINT VARIABLE METHOD 




Since shape of the domain is treated as the design variable, it is 
convenient to think of ft as a continuous medium and utilize the material 
derivative idea from continuum mechanics. The process of deforming H to a new 
domain ft^ may be viewed as a dynamic process as shown in figure 2. One can 
define a transformation as = x + tV(x) where t plays the role of time and x is 
a point in initial domain ft that moves to point x in the deformed domain ft • 

Note that the "shape design velocity" V(x) of point x can be considered as a 
perturbation of design variable. A detailed discussion of this method can be 
found in references 1 and 2. 


The adjoint variable method of design sensitivity analysis (refs. 1, 2, and 
3^ is applied by defining an adjoint equation for an adjoint displacement field 
X to obtain the variation where ft^ is the small region or the finite element 
considered, ^ is a characteristic function for the corresponding ft^, and <J> is the 
normalized von Mises yield stress. 


Note that only boundary integrals appear in the expression for The 
normal movement (V^n) plays the role of shape design perturbation and can be 
expressed in terms of shape design parameters. 



Material derivative of cost: 


J r (V T n)dr 


Material derivative of constraint Adjoint equation 

a(\ k ,» - JJ l [-2*r-(«)o* 3 a)]^ da 

i,j=l do 3 

2 

<!>.=- J r i[ l o 1 -^ (z)e 1 '^(X k ) ] (V T n)dr - J m (ip ~ <|>)(V T n)dr 
k * i,j=l k 


Figure 2 
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PARAMETRIZATION OF BOUNDARY T 


In order to compute design sensitivity the variable boundary should be 
parameterized in terms of a design variable vector b (refs. 2 and 3). Presume 
that points on the boundary F are specified by a vector x (a;b) from the origin 
of the coordinate system to the point S on the boundary, as shown in figure 3, 
where ot is a parameter vector. 

When the vector b of design variables, b = [bp...,b ] , has been defined, 
the domain optimization problem reduces to selection of tne finite dimensional 
vector b to minimize a cost function, subject to the constraints. By defining 
x^(a;b) = x(a;b + TSb), one can define the velocity field at the boundary by 
taking the derivative of x with respect to t. Taking the scalar product of V 
with the unit outward normal to the boundary V and substituting the result into 
the analytical expressions for and ^ yields numerically computable 
sensitivity formulas. 



V [x(a;b + x<5b)] = fib 

(V T n) = [n T — ]6b 

♦o - [/r " T Is dr I 6b 

♦i " l/ r G(z.X k ) (n T ^ )dr] «b 

Figure 3 
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COMPARISON OF FINITE ELEMENT METHODS 

Since shape design sensitivity information is given as a boundary integral, 
one has to check the accuracy of the numerical analysis results on the boundary. 
For comparison of accuracy, constant stress triangular (CST) , linear stress 
triangular (LST) , and 8-noded isoparametric (ISP) elements with optimal stress 
(refs. 4 and 5) are used to calculate design sensitivity. That is, stress values 
are evaluated at Gauss points and linearly extrapolated to obtain boundary 
stresses and strains. 


For boundary parameterization, piecewise linear and cubic spline 
representations are used. In order to compare accuracy of results obtained with 

different finite elements, the same small region should be used to average 

stress. The small regions selected are shown in figure 4, located next to the 

variable boundary where it is most difficult to obtain accurate design 
sensitivity results (ref. 2). 


Define A\k = ^(b + 5b) - (b). The ratio of and Aip times 100 is 
as a measure of accuracy; i.e., T00% means that the predicted change ty' is 
exactly the same as actual change. Numerical results with 6b = 0.001b are 
in figure 4. 


used 

shown 




Sensitivity Check ( ip 1 / Aip x 100)% 


¥ 



Piecewise Linear 

Cubic Spline 

Region No. 

CST 

LST 

ISP 

ISP 

5 

67.5 

99.2 

102.8 

102.6 

6 

68.6 

99.2 

101.8 

101.7 

7 

68.3 

99.1 

100.0 

100.4 

8 

70.1 

99.1 

98.4 

97.4 

9 

79.3 

98.3 

105.2 

104.9 

10 

183.6 

87.0 

102.8 

L 

104.1 


Figure 4 
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OPTIMIZATION OF FILLET 


The cubic spline function, which has two continuous derivatives everywhere 
and possesses minimum mean curvature, is employed here to define the moving 
boundary and 8-noded isoparametric finite elements shown in figure 5 are used for 
analysis. The finite element model contains 131 elements, 458 nodal points, and 
846 degrees of freedom. 

Heights of nodes that define the varied boundary are chosen as the design 
variables, as shown in figure 5. The fillet is optimized using the 
Linearization Method (ref. 6). Convergence criteria require the L-2 norm of 
direction vector p to be zero at the optimum point, where p is obtained by solving 
a quadratic programming problem. For numerical data. Young’s modulus, Poisson’s 
ratio, and allowable yield stresses are 30 x 10^ psi, 0.293, and 120 psi 
respectively. 

The initial design is b = [5.55, 5.1, 4.65, 4.2, 3.75, 3.3, 2.85, 2.4. 

1.95]. _Jnitially, cost, maximum stress violation, and ||p|| are 145.1 in 2 , 

2.1 x 10 , and 2.0_jespectively. A|ter optimization, they are reduced to 

133.4 in , 6.0 x 10 , and 8.8 x 10 , respectively. The final design is shown 

in figure 5 with design variable b « [2.64, 2.13, 1.90, 1.74, 1.61, 1.55, 1.5, 

1.5, 1.5]. 




(a) initial design 


<J> n = 145.1 in. 2 ♦ 133.4 in. 2 


(b) final design 


max ik = 2.1 x 10 ^ 6.0 x 10 ^ 

k 


Figure 5 
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DESIGN OF AN ENGINE CONNECTING ROD 

An engine connecting rod connects the crankshaft and piston pin of an 
engine, transmitting axial compressive load during firing and axial tensile load 
during the suction cycle of the exhaust stroke. The geometry of the connecting 
rod considered is shown in figure 6. Considering that the loads acting on the 
rod are in a plane and that the rod is nearly symmetric about this plane, one 
can reasonably assume that the rod is in a plane stress state. With the main 
interest in the shank and neck regions, the shape of the shank and neck regions 
of the rod are to be determined through the optimization process. The optimum 
thickness distribution, which varies independently from the domain variation, is 
to be determined in the optimization process. To satisfy the condition that the 
distance between the piston pin and the crankshaft is prescribed, it is required 
that the length of the rod not be changed. 



283. S 


Design variables: Thickness h and shape of shank area 

Figure 6 
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DESIGN SENSITIVITY ANALYSIS OF CONNECTING ROD 


The optimal design problem is to find a boundary shape and shank thickness 
to minimize total volume of the rod, with stress constraints ♦ For stress con- 
straints, lower and upper bounds are imposed on averaged principal stresses of 
inertia and firing loads. 

As in the fillet design problem, one can use the principle of virtual work 
to derive a variational equation of elasticity. One can then employ the material 
derivative idea from continuum mechanics and an adjoint variable technique to 
calculate the shape design sensitivity formulas (ref. 7). The sensitivity expres- 
sion resulting from thickness variation can also be found using the same adjoint 
variable method (ref. 2). 

To use the sensitivity formulas computationally, the thickness function h is 
selected to be piecewise constant over strips of finite elements that run along 
the shank. Also, a cubic spline function is used to parameterize the boundary. 
(See fig. 7.) 


Cost: i|> Q = JJ Q h da 

Constraint: ^ = / / <J>(a(z))m^ da < 0 , k = l t 2,»..,2NE 

= J r h(V T n)dr + // fl 6h da 
\ = fj n F 1 (z,X k )6h da + f T F 2 (z,X k )(V T n)dr 


Figure 7 


OPTIMIZATION OF CONNECTING ROD 


With the sensitivity coefficient obtained, one can apply the Linearization 
Method (ref. 6) to obtain the optimum shape and thickness distribution. An 8- 
noded isoparametric element is used for analysis. A finite element model 
including 422 elements, 1493 nodal points, and 2983 degrees of freedom is 
employed For numerical data. Young’s modulus and Poisson’s ratio are 

2.07 x 10 MPa and 0.298 respectively. Upper and lower bounds of principal 
stresses of inertia are 136 MPa and -80 MPa, whereas they are 37 MPa and -279 MPa 
for the firing case. 

The manufacturer’s design is taken as an initial design, where the cost 
functional maximum constraint violation, and ||p|| were initially 726050 mm 3 , 

2.7 x 10 , and 5.9, respectively, and two constraints were active or violated 

around th^ neck area_^near section a— ^). After optimization, they are reduced to 
697182 mm , 1.0 x 10 , and 6.5 x 10 , respectively, with 50 stress constraints 

active. The shape of the initial and final designs and several cross sections 
are illustrated in figure 8. 



Section Q-Q Section b"b Section C - C Section d~d 


- 726050 mm + 697 182 . max ib = 2_7 in ^ 1 ft x in 

Figure 8 
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DESIGN OF A BEAM-PLATE-TRUSS 


Figure 9 shows a truss-beam-plate built-up structure in which thin flat 
plates, stiffened by longitudinal and transverse beams, are supported by four 4- 
bar trusses. A uniformly distributed load is applied to the plates. The points 
supported by trusses are at the intersection of two crossing beams nearest the 
free edges of the structure. The plates and beams are assumed to be welded 
together. The design variable in this problem is the combination of plate 
thickness, width and height of the rectangular cross sections of beams, and 
positions of beams. The design problem is to minimize the volume of the built-up 
structure, subject to constraints on displacement, stress, natural frequency, and 
bounds on design variables. 

The state variable for this built-up structure consists of the plate 
displacements, beam displacements and torsion angles, and nodal displacements of 
the trusses, which satisfy kinematic interface conditions (kinematically 
admissible displacement fields). Hamilton’s principle results in a variational 
formulation of the governing structural equilibrium and eigenvalue (free 
vibration) equations. 
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(b) Side View 


BEAM-PLATE-TRUSS 


Design Variables: 

Beam cross-sectional area 
Plate thickness 
Positions of Beams 

Constraints : 

Displacement 

Compliance 

Eigenvalue 

Stress on beams and plates 


Figure 9 



DESIGN SENSITIVITY ANALYSIS OF BEAM-PLATE-TRUSS 

Design sensitivity analysis with respect to conventional design variable and 
shape using material derivative and adjoint variable method may be extended 
directly to the built-up structure problems. For conventional design variation, 
the general sensitivity formula contains contributions from each structural 
component directly. For shape variation, contributions from each component 
appear as integrals over common boundaries, using interface conditions on the 
common boundaries. 

In figure 10, comparison between actual changes and predictions for 
constraints with 5% changes in all conventional design variable are presented. A 
finite element model of 100 plate elements, 80 beam elements, and 16 truss 
elements is used, with 363 degrees of freedom for total structure. For numerical 
data. Young's modulus, Poisson's ratio, and material density are 3.0 x 10 7 psi, 
0.3, and 0.1 lb/in respectively. Results shown in figure 10 indicate that 
sensitivity accuracy is very good for conventional design. 


= l // ± a F,(z,X k )dh dft + l f F (z,X k )(V T n)dr 

i,j a 3 i,j r iJ 2 

SENSITIVITY CHECK FOR CONVENTIONAL DESIGN 


Constraint 

El. 

No. 

^'/A^xlOO 

Constraint 

El. 

No. 

i|/'/ AiJjxIOO 

Displace- 

C 

112.7 


1 

95.1 

ment 




3 

111.3 


1 

108.8 

stress 

5 

109.5 

Stress 

3 

109.7 

on 

12 

109.7 

on 

5 

109.6 

plate 

14 

109.5 

beam 

11 

106.8 

element 

23 

109.1 

element 

13 

110.0 


25 

109.8 


15 

109.2 


35 

115.1 





45 

113.9 

Eigenvalue 


91.3 



Figure 10 


433 



SHAPE DESIGN SENSITIVITY CHECK FOR BEAM-TRUSS -PLATE 


It is well known (ref. 8) that finite element results on interface 
boundaries, where abrupt changes in the boundary conditions occur (interface 
conditions), are far from being satisfactory. Based on this fact, a finer grid is 
used for shape design sensitivity calculations. Only one quarter of the entire 
structure is used for calculation, due to symmetry. A nonconforming 12 degrees- 4 

of-freedom finite element is used for plates. A finite element model of 400 
rectangular plate elements, 80 beam elements, and 4 truss element is used, with 
total of 1281 degrees of freedom. The same numerical data that are used in 
conventional design sensitivity calculations are used. t 

In figure 11, sensitivity accuracy results are given for 5% uniform changes 
in all shape design variables (positions of beams). Results in figure 11 show 
reasonably good agreement between sensitivity predictions ^ and actual changes 
A\|>k for all except some stress constraints on plate elements. That is, the 
sensitivity results for the stress constraints on plate elements adjacent to the 
interface (marked by *) are poor, even with finer grid. 


Constraint 

El. 

No. 

i|>' / AiJjxlOO 

Constraint 

El. 

No. 

<|>'/ Ai^xlOO 

Displace- 

C 

97.5 


1 

100.7 

ment 




19 

103.6 





44 

109.9* 

Compliance 


92.3 


58 

98.7 




Stress 

85 

188.3* 

Eigenvalue 


95.5 

on 

128 

97.8 




plate 

149 

98.7 


1 

99.6 

element 

175 

91.8 

Stress 

21 

99.7 


198 

109.6 

on 

30 

100.0 


237 

449.0* 

beam 

45 

98.3 


259 

138.8* 

element 

55 

103.1 


296 

95.9* 


77 

103.9 


317 

87.8 





358 

105.3 

i 

■ 



400 

111.3 


Figure 11 
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DESIGN OF A SIMPLE BOX BUILT-UP STRUCTURE 


A simple box built-up structure, in which five plane elastic solid plates 
are welded together, attached to a wall is shown in figure 12. A uniformly 
distributed line load is applied on top of the two side plates and the end 
plate* The shape design variable in this problem is the length d, width b, and 
height h of the box. 

As in the beam-plate-truss case, the principle of virtual work results in a 
variational formulation of the governing structural equations. Then, one can use 
t the material derivative idea and an adjoint variable method to obtain the shape 

sensitivity formula. 

In view of beam-plate-truss shape sensitivity results, an equivalent but 
alternate form of shape sensitivity formula is used for this problem. Since 
finite element results are accurate on the domain and not on the boundary, the 
shape sensitivity formula's expressed in terms of domain integral (refs. 1 and 
2). Hence, instead of (V n) , one has terms V and (div V) in the sensitivity 
formula. 



Design Variables: d, b, and h 


Constraints: i)> k = JJ <|»(a(z))m k dft < 0 , k = 1,2, •••,NE 


% - I // , [gU,z) T V + F(X,z)div V]d8 
K i a 1 


Figure 12 
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SHAPE DESIGN SENSITIVITY FOR SIMPLE BOX 


Since shape design variables are given as d, b, and h, one can assume the 
velocity field to be linear on each plate and thus (div V) is constant. An 8- 
noded isoparametric element is used for analysis. A finite element model of 320 
elements, 993 nodes, and 1886 degrees of freedomyis used. For numerical data. 
Young's modulus and Poisson's ratio are 1.0 x 10 psi and 0.316 respectively. 

The dimension of the structure is b = d = h = 8 in. and the thickness of the plates 
is 0.1 in. Uniform external load is 4.77 lb/in. 

In figure 13, the sensitivity accuracy result is given separately for 3% 
change in d and h. Results given in figure 13 show excellent agreement between 
predictions and actual changes A\|»^. The boundary method that is applied to 
the beam-plate-truss built-up structure is tested to the same box problem with 
unacceptable results. The domain method of shape design sensitivity for built-up 
structure has a promising future. Work continues in evaluating the method on 
larger scale examples. 


Region 

El. 

No. 

x 100 

6d = 0.03d 

x 100 

6h = 0.03h 

Region 

El. 

No. 

4? 1 / A\J> x 100 
6d = 0.03d 

Tp' / Alp X 100 
6h = 0.03h 


1 

97.6 

103.6 


129 

96.2 

103.4 


15 

98.2 

103.6 


143 

98.1 

103.0 

Top 

29 

98.8 

103.5 

Side 

157 

98.0 

102.6 


43 

100.4 

103.6 


171 

95.6 

103.0 


57 

110.1 

103.5 


184 

100.4 

103.2 


72 

97.9 

103.3 


264 

99.9 

103.7 


86 

98.5 

103.2 


278 

114.1 

103.9 

Bottom 

100 

98.5 

103.1 

End 

291 

95.1 

101.1 


114 

100.9 

103.3 


307 

120.3 

103.5 


128 

100.9 

103.2 

1 

320 

99.9 

103.7 


Figure 13 
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